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The kinetics of geminate subdiffusion-assisted reactions (SDARs) of interacting particles is an- 
alyzed in detail with the use of the non-Markovian fractional Smoluchowki equation (FSE). It is 
suggested that the interparticle interaction potential is of the shape of potential well and reactivity 
is located within the well. The reaction kinetics is studied in the limit of deep well, in which the FSE 
can be solved analytically. This solution enables one to obtain the kinetics in a simple analytical 
form. The analytical expression shows that the SDAR kinetics fairly substantially depends on the 
mechanism of reactivity within the well. Specific features of the kinetics are thoroughly analyzed in 
two models of reactivity: the subdiffusion assisted activated rate model and the first order reaction 
model. The theory developed is applied to the interpretation of experimental kinetics of photo- 
luminescence decay in amorphous a-Si:H semiconductors governed by geminate recombination of 
electrons and holes which are recently found to undergo subdiffusive relative motion. Analysis of 
results demonstrates that the subdiffusion assisted activated rate mechanism of reaction is closer to 
reality as applied to amorphous a-Si:H semiconductors. Comparison of experimental and theoreti- 
cal kinetics allowed for obtaining some kinetic parameters of the systems under study: the rate of 
escaping from the well and the parameter characterizing the deviation of the subdiffusive motion 
from the conventional one. 
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I. INTRODUCTION 

Diffusion assisted reactions (DARs) is the important 
stage of a large number of chemical and physical con- 
densed phase processes jii^i^ In many of these processes 
the DARs are known to be strongly affected by inter- 
action between reacting particles. The effect of the in- 
teraction on the DAR kinetics is actively studied both 
experimentally and theoretically for a long t\me^^ 

In the case of conventional diffusion the problem of the 
analysis of DAR kinetics reduces to solving the Smolu- 
chowski equation for the pair distribution function (PDF) 
of reacting particles. This equation is fairly complicated 
and can, in general, be solved only numerically. As for 
analytical study, usually it is made with the use of steady 
state analytical solutions.— The general time dependent 
analytical solutions can be found only for very few in- 
teraction potentials, for example, in the case Coulomb 
interactioni^ This solutions, however, are very cumber- 
some and complicated for applications. 

Some years ago much more simple and rigorous method 
of analytical solution of the Smoluchowski equation was 
proposed, which is applicable in the practically interest- 
ing limit of deep well of attractive interaction (reactiv- 
ity assumed to be localized in the well) The solution 
shows that the interaction strongly manifests itself in the 
reaction kinetics resulting in the long life time of parti- 
cles within the well (i.e. caging). The time evolution 
of the PDF of pairs captured and reacting in the well 
appears to be non-exponential.^'^ This specific feature 
of the PDF time evolution shows itself, for example, in 
non-exponential kinetics of geminate DARs with the long 
time tail of inverse power typei^ii 



In this work we will consider the kinetics of gemi- 
nate reaction of interacting particles undergoing subd- 
iffusive motion. Recall that subdiffusion is a certain 
type of anomalous diffusion, which is characterized by 
the anomalously slow time dependence of the mean 
square of displacement {r'^(t)) ~ with a < 
Recently, the specific features of the kinetics of subd- 
iffusion assisted reactions (SDARs) is a subject of ac- 
tive discussions^^'i^^'ii^'i^^i^^^ The anomaly 
of diffusion is shown to affect fairly strongly the reac- 
tion kinetics leading to the effective slowing down of the 
reaction at long times, to the strong fluctuations of con- 
centrations of reacting particles at long times, etc. 

In the absence of interparticle interaction (i.e. in the 
case of free subdiffusion) the time evolution of the PDF of 
subdiffusing particles is usually described by the analog 
of diffusion equationj^i^ which is called fractional diffu- 
sion equation and in which the effect of diffusion anomaly 
shows itself in anomalously long time memory. The pre- 
dictions of the theory based on fractional diffusion equa- 
tion are analyzed in a large number of papers (see, for 
example, reviews [8] and [9]). 

As for the SDARs of interacting particles, these pre- 
cesses are not studied theoretically yet. The kinetics of 
them is also determined by the corresponding PDF, but 
the PDF evolution is described by the fractional Smolu- 
chowski equationj-^'^^ Similarly to the case of conven- 
tional diffusion, the factional equation is much more com- 
plicated for numerical and analytical analysis than that 
for free diffusion. In this work we propose the analytical 
solution of the fractional Smoluchowski equation in the 
limit of deep well, assuming the reactivity to be localized 
within the well. With the use of the obtained solution 
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the PDF evolution and the kinetics of geminate SDARs 
are analyzed in detail in this limit. 

The analysis shows that, unlike the case of conven- 
tional DAR, the SDAR kinetics strongly depends on the 
mechanism of reaction. In particular, the kinetics ap- 
pears to be essentially different for two models of reac- 
tivity: the subdiffusion assisted activated rate model [or, 
more generally, the kinetically (i.e. mobility) controlled 
reaction model] and the first order reaction model. This 
strong difference enables one to select the realistic reac- 
tion mechanism by comparison of theoretical predictions 
with the experimental data. 

The obtained results are applied to the interpretation 
of the experimental kinetics of photoluminescence decay 
in amorphous semiconductors a-Si:H resulting from gem- 
inate recombination of photoexcited electrons (e) and 
holes (/i)i^i2ii2^ Electrons in these semiconductors are 
known to be highly mobile, undergoing subdiffusive (dis- 
persive) migration, while holes are nearly immobile!^ 
Recently, fairly detailed experimental investigation of the 
kinetics of geminate e-h recombination at different tem- 
peratures has been carried out^^ and experimental re- 
sults have been semiquantitatively described within the 
free subdiffusion model. It is worth noting, however, that 
the kinetics of the process under study is, clearly, signif- 
icantly affected by the Coulomb e-h interaction which is 
quite strong in the investigated semiconductors: for di- 
electric constant e « 10 characteristic for these semicon- 
ductors and temperatures T < 300K the Onsager radius 
(distance, at which the Coulomb interaction is equal to 
the thermal energy) is estimated as le > 50A. 

Analysis of theoretical kinetic dependences, obtained 
in this work, shows that the activated rate model de- 
scribes the experimental results better than the first or- 
der reaction one. This analysis allowed us to obtain the 
characteristic model parameters of the system: the pa- 
rameter a characterizing the diffusion anomaly, the rate 
of escaping from the e-h interaction well, etc., which re- 
sult in the best fitting, and estimated their dependence 
on temperature. The proposed interpretation is com- 
pared with another one applied in ref. [24] to treat some 
earlier experimental results on the same process. 



II. DIFFUSION ASSISTED REACTIONS 

We start our analysis with the discussion of gemi- 
nate reactions of interacting particles assisted by con- 
ventional diffusion. The reacting particles are assumed 
to be spherically symmetric. We also assume that the in- 
terparticle interaction potential is spherically symmetric: 
u{r) = u{r) = U{r)/kBT, where r = |r| is the interpar- 
ticle distance, and is of the type of potential well with 
the reaction barrier at short (contact) distance r = d, 
the bottom at r — ri, (see Fig. 1). It is suggested, in 
addition, that the well is deep enough so that the ac- 
tivation energies of escaping from the well and reaction 
are large: Ua = —u{rb) ^ 1 and u{d) — u{rb) ^ 1. In 




FIG. 1: Schematic picture of the distance dependence of the 
interaction potential u{r) = U {r) / (ksT) and reactivity kr{r). 
In this picture d is the distance of closest approach, rj, is the 
coordinate of the bottom, and U is the Onsager radius. The 
dashed line shows the Coulomb potential which leads to the 
cusp- like well at r = d. 



this limit the characteristic time of reaction in the well 
is much larger than the time of equilibration in the well 
and the long distance (Coulomb) part of the potential 
can be characterized by the effective Onsager radius 



r°° dr r 



'2gU(r) 



(2.1) 



nearly independent of [this radius satisfies the relation 
u{le) ~ 1]. Noteworthy is that in the deep well limit the 
explicit shape of the well of u{r) at short distance turns 
out to be not important for the kinetics: it can be, for 
example, of parabolic shape with ri, > d or edge type one 
with Tb d (see Fig. 1). 

In the Markovian approach implied by the conventional 
diffusion approximations the DAR kinetics is described 
by the PDF p(r, t) of reactive pairs, satisfying the kinetic 
equation 



p{v,t) = -Cp{r,t), 



(2.2) 



where £ is the operator which determines space/time 
evolution of the system under study. The form of this op- 
erator depends on the process considered (see Sec. III). 

In this work the observable under study is the geminate 
DAR kinetics, i.e. the time dependent probability Yr{t) 
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of geminate reaction, which we will call the DAR yield. In 
our further analysis, however, it will be more convenient 
to analyze the DAR flux 

Jr{t)=Yr{t). (2.3) 

To obtain the expression for DAR yield and DAR flux 
one needs to specify the model of reactivity. In our work 
we assume that the reactivity is localized within the well 
and consider two most well known models of reactivity: 
the diffusion assisted activated rate and first order reac- 
tion models. They correspond to two different reaction 
mechanisms: kinetically controlled (controlled by relative 
migration) and first order reaction controlled reactivity, 
respectively. Possible examples of these reaction models 
are discussed below as applied to analysis of experimen- 
tal results (see Sec. IV). Here we will restrict ourselves 
to discussion of the mathematical details of the models. 

1) Activated rate model. This model treats the DAR 
as a diffusive flux over a barrier located at the reaction 
surface [in the studied model of spherically symmetric 
particles this is the barrier at r = d (see Fig. 1)]. Some 
well known formulas for the reaction yield obtained in 
this model are presented below (in Sec. IIB). The ac- 
tivated rate model is, actually, a particular example of 
the general class of kinetically controlled reaction mod- 
els, which predict the reaction rate proportional to the 
diffusion coefficient, or more generally speaking, to the 
mobility of particles (Sec. IV). 

2) First order reaction model. In this model the DAR 
flux Jr{t) is determined by the fluctuating rate kr[r{t)\ 
of first order reaction and is expressed in terms of the 
functional: 

Jr{t) = (krW)] exp ^ dTkMr)]j^ , (2.4) 

in which the average is made over the stochastic fluctua- 
tions of r{t) governed by process (|2.2[) . 

In the considered approach the DAR flux can be writ- 
ten in a universal form in terms of the the Green's func- 
tion p{r,ri\t) of the stochastic Liouville equation 

p = -[C + kr{r)]p + 6{t)5{r - r,) : (2.5) 
Mt) = - [ dv [ dr,p{v,r,\t)p,{Yi), (2.6) 

Jr>d J 

where Pi{r) is the initial spatial distribution of particles 
[naturally, in the activated rate model kr (r) = 0] . 

A. Mechanisms of stochastic motion 

Here we will briefly discuss some useful models for de- 
scription of relative jump-like stochastic motion of re- 
acting particles and analyze the validity of the diffusion 
approximation for description of the process under study. 



One of the most general models of spatial jump-like 
evolution of the system is based on the Kolmogorov- Feller 
approach in which 

C^~w(l~V), (2.7) 

where w is the jump rate, for simplicity assumed to be in- 
dependent of the coordinate, and P is the matrix of tran- 
sition probabilities satisfying the detail balance relation 
and the normalization condition which in the coordinate 
representation V{Y,Yi) for V is written as 

/dr7'(r,r0 = 1. (2.8) 

In general, there is a large variety of jump models for 
V satisfying above relations. Here we will discuss a class 
of models especially convenient for theoretical analysis. 
These models are based on the representation of the ma- 
trix V in terms of the differential Smoluchowski-like dif- 
fusion operator 

Lp = -I?Vr(Vrp-|-pVrU), (2.9) 

in which Vr is the gradient operator in {r}-space and D 
is the diffusion coefficient f^^i^l 

/•OC 

V^P{L)^ dre-^^Pir). (2.10) 
Jo 

In the representation (|2.10p the function P(t) is properly 
normalized: dr P{t) = 1, so that V can be consid- 
ered as the operator of diffusive evolution averaged over 
some distribution function P{t) of evolution times. Evi- 
dently, the operator V (|2.10p , which has the form of the 
Laplace transform of -P(t), satisfies the normalization re- 
lation ((2?8| . 

The model (|2.10p is very useful for the analysis of ap- 
plicability of the diffusion (Smoluchowski) approximation 
which appears to be valid in a wide region of parameters 
of the model and times. The validity criterion can eas- 
ily be obtained with the use of eq. (|2.2p written for the 
Laplace transform in time p(r, e) = dr p{r,T)e~'^'^ : 

ep = -w{l-'P)p. (2.11) 

According to eq. (|2.1ip for small e/w <C 1, when the left 
hand side of this equation is small, the operator I — V in 
the right hand side can be expanded in Lt: 1 — V~ —Lt, 
where 

i = dttP{t) (2.12) 

is the average time of diffusion-like evolution resulting in 
the jump operator V. The effect of next orders of expan- 
sion of V in Lt is evidently of higher orders in e/w <C 1 
and therefore is negligibly small. The correctness of this 
statement can also be demonstrated by expanding the 
solution of eq. (|2.1ip in the (complete) basis of eigen- 
functions of the operator L. 
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The presented analysis shows that at relatively long 
times t > 1/w the kinetics of processes, governed by 
jump-like operator (|2.10p . is quite accurately described 
by the corresponding Smoluchowski equation. 

In what follows we will restrict ourselves to this long 
time limit of the reaction kinetics and correspondingly to 
the Smoluchowski approximation. In addition, for sim- 
plicity, we will consider the spherically symmetric gemi- 
nate processes which are described by distribution func- 
tions depending only on distance r — |r|. 

B. Equations of diffusion approximation and 
two-state model 

In the considered case of spherically symmetric gemi- 
nate reactions the PDF p(r, i), depending on the inter- 
particle distance r — |r|, satisfies the Smoluchowki ki- 
netic equation, which for 

a[r,t)^rp{r,t) (2.13) 

is written as 

CT = £)Vr(VrCr ctV^m) - fc^cr, (2.14) 

where = d/dr and D is the diffusion coefficient. 

The function a{r,t) (|2.14[) obeys the boundary condi- 
tions a{r ^ oo,t) — + and a{r — > 0,t) 0, and the 
isotropic initial condition 

a,{r) ^ a{r,0) = S{r - r,)/{4TTn) (2.15) 

with Vi ri, corresponding to the creation of parti- 
cles within the well. This equation can also be rep- 
resented in terms of the Laplace transform cr(r, e) — 

ea — Ui = _DVr(V.rCr -|- crVfU) — kr^. (2-16) 

This representation appears to be more convenient for 
our further analysis. 

In general, eqs. (|2.14p and (|2.16p cannot be solved an- 
alytically. In the limit of deep well, however, the solution 
can be obtained in a simple form by expansion in a small 
parameter r^/re <^ 1, where ~ ^t/D is the time of 
equilibration within the well and Te ^ r^e^"" is the time 
of escaping from the wellj^ii 

Analysis of this solution shows^ that in the lowest or- 
der in the parameter r^/Te the Smoluchowski approxima- 
tion (|2.14p and (|2.16p is equivalent to the model of two 
kinetically coupled states. This two-state model treats 
the process under study as transitions between the state 
within the well (d < r < l^), whose population is 

n(i)=47r/" drra{r,t), (2.17) 
Jd 

and the state of free diffusion outside the well {r > I) de- 
scribed by the distribution function C{r,t) = r~'^c{r,t). 



The functions n{t) and c(r, t) satisfy simple kinetic 
equations,^''' which can be written in the most compact 
form in terms of Laplace transforms n{e) and c(r, e): 

en ^ 1 + [S^^K+c{le, e) - (K^ + Wr)n] (2.18a) 
ec = [DVlc+ {SiK^n- K+Z)d{r (2.18b) 

where Si — {AttI^)^^. The terms proportional to K± de- 
scribe the above-mentioned kinetic coupling (transitions) 
between the state within the well, located at r = l^, and 
the free diffusion state outside the well. In the consid- 
ered limit Tr/Te <ti 1 the transition rates K± satisfy the 
relations?^ 

K±^oo and K+/K_ =K^ = (2.19) 

where 

Z^^ f drr^e-"^'') (2.20) 

is the partition function for the well. 

Equations (|2.18p are written for the initial condition 

n(0) = 1 and c(r, 0) = 0. (2.21) 

corresponding to the initial population of the well implied 
by eq. (|2.15p . As to the boundary conditions for c(r, i), 
they are given by ZeVrc(r, t) — c{r, t)\r=i^ = and c(r 

oo) = 0. 

The term WrU in eq. (|2.18ap describes the effect of the 
first order reaction in the well. In the considered (Marko- 
vian) diffusion approximation the two above-mentioned 
models of reactivity in the well result in the similar ki- 
netic equations of the form (|2.18p . The only difference 
consists in the analytical expression for Wr-: 

a. Activated reaction model (kr = 0). In the dif- 
fusion assisted activated reaction model, implying acti- 
vated diffusive passing over the barrier at r ~ d, one 
gets^i 

w^=Wr^ = ^( f drr-2e"(''A . (2.22) 

As we have already mentioned above this model is an 
example of a large class of kinetically controlled reaction 
models, in which the reaction rate is determined by the 
mobility of particles: Wr ^ D. Some examples of such 
models are discussed in Sec. IV. 

b. First order reaction model {kr 0) ■ In the model 
describing the reaction as a first order process with rate 
kr(r) the expression for , for example, in quite realistic 
case of relatively small values of kr within the well^:-'^ 

Wr = Wrj = (kr) = — [ dr kr{r)r^ e'"'^'''^. (2.23) 

Jw 

In this equation the parameter w, used as a limit of in- 
tegration, denotes integration over the region near the 
bottom of the well {d < r < l^)- 
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In both models the Markovian DAR kinetics, expressed 
in terms of the Laplace transform Jr{e) of the DAR flux 
p.3|) . is proportional to the the Laplace transform n{e) 
of the well population: 



Jr(e) = Wrn{e), 



(2.24) 



i.e to analyze the reaction kinetics one should find the 
time dependent well population n{t). 

For our further discussion it is convenient to represent 
equations (|2.18p in a matrix form: 

eR= -(A + W'^)R + R„ where R = (n,c)^, (2.25) 

Ri = (IjO)^ is the vector representation of the initial 
condition ^2^, 



A = 
and 



K_ -Si^^K+JdrS{r-le) 



Wr 



Wr 





(2.26) 



(2.27) 



This representation shows that equations (|2.18p can 
be considered as a two-state analog of the SLE (|2.5p . 
corresponding to the two-state model of fc^- fluctuations 
in the functional (|2.6[) . Naturally, in this SLE the effect 
of reactivity is represented by (2 x 2)- matrix of the form 

It is also worth emphasizing that the two-state model 
(I2.18P is actually based on the approximate replacement 
of the Smoluchowski operator Lr — I?Vr(Vr + V ru) by 
the simpler one A + Wr, operating on the reduced PDF 
R. The approximation is valid in the limit of deep well, 
when all eigenvalues of representing population re- 
laxation within the well arc much larger than the low- 
est one, which describes quasistationary escaping from 
the well and reaction. In the deep well limit the effect 
of high eigenvalues is negligibly small and the kinetics 
of the process is quite accurately treated within the ap- 
proximation taking into account the coupling of the low- 
est state in the well with the continuum of states outside 
the well (which is just equivalent to the proposed two- 
state moder). This main idea of the proposed method 
is important point in our further analysis of SDARs. 

It is also important to note that in the deep well limit 
the two-state model is valid independently of the shape 
of the well near the bottom of the well (see Fig. 1). 
The shape manifests itself only in the value of partition 
function defined in eq. (|2.20p . 



C. Reaction kinetics 



In accordance with the relation (|2.24p the DAR kinet- 
ics [i.e. the reaction flux Jr{i)\ is determined by that 
of the well depopulation n{t) which can be obtained by 



solution of eqs. (|2.18p and subsequent inverse Laplace 
transformation:^''' 



^oo+0 

, , 1 f exp(et) 

— -ioc+O 



(2.28) 



The function V{e) is directly related to the Green's 
hmction of the operator which controls diffusion outside 
the well [with the reflective boundary condition {le^g — 



g{r,r,\e)~{r\{e + DWl) |r,) : 



(2.29) 



Vie) := l/g{k, Qe) = D[l-' + {e/DY'\ (2.30) 



Substitution of the expression ()2.30|) into eq. ([22 



leads to the following formula for the well population 



n{t) 



1 

27ri 



■'°°+"^^_exp[£(^„t)] 



— 'tOC-|- 



e+*fKi)-e-*iK*) 



(2.31) 



where 7 = y/{we/w„){l'^We/D), e±^ i\Jl- jl'^ 



are the roots of equation z 



7z - 



1 = 0, and 



$f (z) = [1 - erf(e±V^)] exp(e±2:). 



The rate 



We 



Wr 



(2.32) 



(2.33) 



is a sum of the rate of escaping from the well 



D 



00 



DIJZ^ (2.34) 



and the rate of reaction in the well Wr [see eqs. (j2.22p 
and ([2:231) ]. 

Specific features of the well population n{t) (|2.3ip are 
analyzed in detail in refs. [6] and [7]. In general, this 
function is non-exponential and the analytical properties 
are essentially determined by the parameter 7 introduced 
in eq. ()2.3ip . The physical meaning of this parameter is 
clear from the relation 7 = -y/ {we/wg){leWe/D) ~ h/lo, 
in which I £, = ^/D/wo is the average distance of diffusive 
motion during the life time tq — Wq of the particle in the 
well. This relation shows that if the well is deep enough. 
Id ^ le and therefore 7 <C 1. The parameter 7 controls 
the qualitative change of the analytical behavior of the 
kinetics: 

n{t) = e^p{-w„t) atT<ln(l/7), (2.35) 
n{t) - 1/^3/2 atr>ln(l/7). (2.36) 
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In the limit of small 7 <C 1 more detailed analysis of 
the dependence n{t) is possiblej^ii however, here we will 
restrict ourselves to these simple relations only. 

In addition to the kinetic time dependences the steady 
state characteristics of type of the total DAR yield Y^r^ — 
Yr{t — > 00) are also of certain interest for further appli- 
cations. They can easily be obtained with the use of 
expressions derived. For example, 



dtJr{t) = 'Wrn{0) = — 

Wr 



■ We 



(2.37) 



III. SUBDIFFUSION ASSISTED REACTIONS 

A. Equations in subdiffusion approximation 

In the case of subdiffusive motion the evolution of 
p(r, t) is described the subdiffusion (fractional) variant 
of the Smoluchowski equationi^ This equation can be de- 
rived within the continuous time random walk approach^ 
for jump-like motion of particles assuming long time 
tailed behavior of the probability density function of 
waiting times for jumps and using formula (j2.10p for 
the distribution of jumps P{r,ri). Similar to the case of 
the conventional diffusion considered (Sec. II A) the frac- 
tional Smoluchowski equation is obtained in the limit of 
relatively weak deviation of p(r, t) from the equilibrium 
PDF, when {tL)p ^ p and therefore one can expand 

V = L in powers of Lt: 



(3.1) 



where = d/dr , Dg is the subdiffusion coefficient, and 

is the Riemann-Liouville fractional derivatives^ with a < 
1. The function a{r,t) (|3.1[) satisfies the boundary and 
initial conditions similar to those used for the case of 
conventional diffusion [see Sec. IIB]. 

The corresponding equations for the Laplace transform 
a(r, e) is written as 

ea - C7, = Dse^^^VrC^r^ + ^'^ru). (3.3) 

It is seen that from mathematic point of view in the 
absence of reaction the difference of the subdiffusion vari- 
ant of the Smoluchowski equation for Laplace transform 
a from the conventional Markovian equation reduces to 
the replacement of D with ZJ^e^^". The problem be- 
comes more complicated, however, when one is going to 
analyze the effect of reactivity. Below we will obtain the 
kinetic equations describing the SDAR processes in the 
two models of reactivity proposed above: the activated 
rate model and the model of first order reaction. 

Since the coordinate part of the subdiffusion differen- 
tial equations (|3.ip and (|3.3p coincides with that of the 



Smoluchowski equation (I2.14p and (|2.16p the two state 
model is still applicable though with some modification 
of the form of time dependences or, as applied to equa- 
tions for Laplace transforms, the dependences on the pa- 
rameter e [i.e. the functions n(e) and c(r, e)]. 



B. Activated rate reaction in the well 

In the model of subdiffusion assisted activated rate re- 
action in the well the corresponding fractional analog of 
the kinetic equation (|2.25p for R = (n,c)^ can be writ- 
ten by taking into account that the fractional Smolu- 
chowski equations (|3.2p and (|3.3p differ from the conven- 
tional ones only in operators which determine the time 
dependence of the PDF. As for the coordinate opera- 
tors, they are similar in both equations differing only 
in diffusion coefhcients: Lr^ = I?sVr(Vr -I- V^-u) = 
DV riy r + ru) ■ This fact is important because the two- 
state model actually reduces to the special two-state rep- 
resentation of the coordinate operator. Correspondingly, 
if in the deep well limit this model is valid in the case 
of conventional diffusion it is also valid as applied to the 
subdiffusion processes since the validity conditions (em- 
phasized in the end of Sec. IIB) appeal only to the char- 
acteristic properties of the coordinate operator (eigenval- 
ues, describing population relaxation in the well, should 
be mach larger than the lowest eigenvalue representing 
the escape rate). The similarity of coordinate operators: 
Lr^ = {Ds/D)Lr, enables us to easily obtain the corre- 
sponding subdiffusion variant of the two-state represen- 
tation ofLr/. As + W^^ = {Ds/D){A + Wr), in which we 
set Wr^ = [{D s / D)Wr]^ ^ to reveal the anomalous di- 
mensionality of Ds- Thus, finally, we arrive at equation 

[e -I- ei-"(A, + W^J]R = R^, where R = (n, 2)^, (3.4) 
Ri = (1,0)^ is the initial condition [see eq. (|2.25p ]. 



A, 
and 



-Sf^K^JdrS{r-Q 



Wr, 





(3.5) 



(3.6) 



is the matrix of SDAR rates in the well. 

The parameters in these equations are similar to those 
in eqs. (^1^ [see eqs. (PTg|) and ^^], except the effec- 
tive rate Wr, of (anomalous) reaction in the well modeled 
by subdiffusive passing over the barrier at r ~ d: 



Wr, 



£1 



l/a 



(3.7) 



Recall that in our calculation we take the limit (j2.19p . 
i.e. K^-^oo with Ke = K^/Kt = Z^. 



7 



Solution of eqs. (|3.4p and inverse Laplace transforma- 
tion give for the well population 



1 



ioo+0 



de- 



exp(et) 



— ioc+O 



(3.8) 



As in the case of conventional diffusion the function 
V{e) is determined by the Green's function of the oper- 
ator describing subdiffusive relative motion of particles 
outside the well [with the reflective boundary condition 

{leVg- g)\r=i, =0] 

gir,n\e)^{r\{e'' + D,Viy'\n) : (3.9) 



V{e) = l/5(/e, le\e) = DS;' + (e7^s)^/']. (3.10) 



Formulas (|3.8p - p.lOp reduce the problem of calculat 
ing n(t) to the evaluation of the Green's function 
and subsequent inverse Laplace transformation ([3^ 



nit) 



ioo+O 



ioo+0 



1 

2TTi 



where 



(3.11) 



(3.12) 



the parameters = ± 1 — jJa ^-''^ similar to 
those introduced earlier [see eq. 

^M], and $±(r) 
are expressed in terms of the Mittag-Lcffler function 



ioo+0 



exp[£(w„^t)] 



de 



(3.13) 



It is easily seen that for a = 1 $^ (t) coincide with func- 
tions defined in eq. (|2.32p . Similarly to formulas for 
conventional diffusion in eq. (j3.1ip we have introduced 
the rate 



I/O 



in which 



(3.14) 



(3.15) 



is the rate of subdiffusive escaping from the well and Wr, 
is the anomalous SDAR rate in the well [see eqs. p.7p ]. 

Unlike the conventional diffusive well depopulation 
kinetics (|2.3ip . the subdiffusive kinetics p. lip is non- 
exponential at all times. Moreover, in the most inter- 
esting case of small 7 <C 1 with high accuracy one can 



neglect the last term in denominator of the integrand in 
eq. p. lip and get the expression 



n{t) ^ Le,[-{w„tr], 



(3.16) 



which predicts the behavior n{t) ~ {wgt)'" at Wgt ^ 1. 

With the use of relation (|2.24p . which in the subdiffu- 
sion case is written as 



Jr{e) 

formula similar to eq. 
SDAR flux Jrit) 



^<n(e), 



(3.17) 



can also be obtained for the 



j(t) = 



exp[£(w, t)] 
de 



ioo+O 

(wr„/woJ" d 



a/2 



dt 



e-_<^>l{t) ~ el^-{t) . (3.18) 



This expression shows the important peculiarity of the 
activated variant of the SDAR kinetics: the time de- 
pendence of the SDAR flux [t) is essentially different 
from that of the well population n{t). Unlike relations 
obtained above for the conventional diffusion and pre- 
diction of the expression p.l6p . in the considered limit 
7 ^ 1 at relatively short times r < 1/7^/" the flux 
Jr{t) ~ La[—{w[)j[,)°'] so that the time dependence of the 
SDAR flux is of inverse power type: 

Jr{t) - {wot)-^^+°'^ for 1 < Wat < 7"^/", (3.19) 
(<) - (u;o/)"^'+"/'^ foru;o/>7"'/". (3-20) 



J, 



Interestingly the expressions for steady state charac- 
teristics like the total SDAR yield = Yr, (t 00) 
are fairly similar to those derived above (see Sec. II) for 
the case of conventional DAR, though with Wr and We 
replaced by w"^ and w"^ , respectively. In particular, 

r°° ~ w" 

r-=/ dtMt)^M0)^ a 1: a - (3-21) 



C. First order reaction in the well 

In the case of first order reaction in the well the 
SDAR kinetics is described by more complicated non- 
Markovian stochastic Liouville equation^i^ii which is the 
non-Mar kovian analog of eq. (|2.14p . In what follows, 
however, for simplicity we will restrict ourselves to the 
more simple two-state model whose subdiffusive variant 
is represented by eqs. p.4p . With the use of the method 
recently proposed in ref. [7] the corresponding two state 
non-Markovian kinetic (stochastic Liouville) equation for 

the vector R = (n,c)^ can be derived in the following 
form^ii^s 



{n + Asn^-°')Ti^R„ where R=(n,c)" 



(3.22) 
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Ri = (1, 0)^ is the initial condition, 
Cl = e + Wrj , and Wr^ — 



Wrf 





(3.23) 



with Wrf defined in eq. (|2.23p . 

Solution of this equation (in the limit K'^ — > oo with 
K"/K" = Z^) yields for the population of the well 



nit) 



where 



2m 



l + £"+7^.(e-e^)"/2' 



w, 



and 7, = {llw'^jD.Y/^. 



(3.24) 



(3.25) 



Naturally, in the model of first order reaction the SDAR 
flux Jr{t) is proportional to the well population: 



Jr{t) — WrfTlit) 



(3.26) 



similarly to the case of conventional diffusion. 

Expressions p.24p and p.26p show that in the first or- 
der reaction model the SDAR kinetics differs from that 
obtained in the activated rate model [eqs (|3.11|) and 
(|3.18p ]. In particular, in the most interesting limit of 
small 7, ^ 1 and for weak reactivity Wr ^ We^ at short 

times t < l/ijl^^WeJ 



J^.{t) = Wrn{t) W Wr e '^''f*La[-{Wej)"] 



(3.27) 



More detailed analysis allows one to conclude that for 
7^ <C 1 over a wide region of times the time dependence 
of the SDAR flux is of inverse power type given by: 

Jrit) ^ t^" iOT 1 < Wet < "f~'^^°' ,WeJWrf\ 

Jr{t) - i-(l+"/2) for ^ WeJWrj, ij^'""- 

In other words, the long time behavior of Jr{t) is simi- 
lar in both models of reactivity, however at intermediate 
times the first order reaction model predicts slower ki- 
netics than the activated rate one. 

It is of certain interest to note that despite this dif- 
ference in the SDAR kinetics the formula for the SDAR 
yield in the first order reaction model appears to coincide 
with that in the the activated rate model [see eq. p.2ip ]: 

i;°° = <^/«^-hO. 



IV. ELECTRON-HOLE RECOMBINATION 
KINETICS 

Here we will apply the obtained results to describing 
the kinetics of photoluminescence decay in amorphous 
semiconductors a-Si:H governed by subdiffusion assisted 
geminate recombination of photogenerated charge car- 
ries: electrons (e) and holes {h)i^-2^^ The geminated 



e-h recombination is strongly influenced by the attrac- 
tive interaction u(r) of e- and /i-quasiparticles. At e-h 
distances r much larger than the size leh of these quasi- 
particles the potential u{r) is of the Coulomb form 



u(r) 



-le/r, with le ^ e^/^esksT), 



(4.1) 



where le is the Onsager radius [see eq. (|2.ip ] and eg « 10 
is the static dielectric constant for the semiconductor a- 
Si:H.^'^'^^ Note that leh is probably of order of the mean 
spacing in this semiconductor. Is is also worth noting 
that at short distances r ~ rf, <C leh the potential is ex- 
pected to be somewhat flattened because of the effect of 
finite size of charge distribution at these distances result- 
ing in the finiteness of e-/i-interaction at r <C leh- This 
effect can qualitatively be understood by considering the 
potential of interaction between point-like charge with 
the homogeneously changed sphere which is known to be 
of parabolic shape and finite at the center of the sphere. 
In the case of quantum particles the qualitatively same 
effect is also expected as one can see form the analysis 
of the electronic terms of the simplest molecules H2 and 

However, independently of behavior of the potential at 
short distances this potential is of the shape of the po- 
tential well (see Fig. 1) and in the deep well limit the 
theory proposed is valid independently of the well shape, 
as it has been mentioned above. The characteristic fea- 
tures of well shape manifest themselves only in the value 
of the partition function Z^u and thus in values of rates 
Wr and We which are considered as adjustable parameters 
anyway. 

It is worth noting that the Onsager radius le is fairly 
large for the systems considered (with the dielectric 
constant e « 10) for the temperatures of experiments, 
100 < T < 170 (K), we get 140 > le > 80(A). In such 
a case the assumption that the geminate e-h pairs are 
initially created within the well, i.e. sd r < le looks 
quite reasonable. The large value of the Onsager radius 
le ^ i"h ^ leh implies the large depth of the potential 
well. In addition, large value of le also ensures the valid- 
ity of the diffusion approximation for description of the 
jump like motion of electrons applied in our analysis. 

The exact shape of the potential and the mechanisms 
of jump-like motion at short distances r ^ d < leh, which, 
according to the above discussion, determines the mecha- 
nism of reaction in the well, can hardly been found explic- 
itly. One of the goals of this section is the selection of the 
proper model of recombination process in the well, i.e., 
actually, selection of mechanism of reactivity, by com- 
parison of predictions of the models with experimental 
results. In our analysis we will consider two of them dis- 
cussed above: 

a) First order reaction model. This model assumes that 
e-h recombination in the well is the first order reaction 
whose rate is determined by the direct charge-transfer 
exchange interaction^ [the model implies the exponen- 
tial distance dependence of the reaction rate kr{r) ~ 
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exp(— r/re)]. This model is traditionally applied to de- 
scribing the kinetics of condensed phase recombinationi^ 

b) Activated rate model. The activated rate model is 
based on the mechanism of diffusive passing over a barrier 
which can be considered as kinetically controlled reaction 
(i.e. reaction controlled by mobility of particles). This 
model is quite realistic as applied to charge recombina- 
tion in polar liquids in which the mean-force interaction 
potential of charges is known to have a fairly high po- 
tential barrier at short distances r ~ d.^^ As for solids, 
one can hardly expect similar barrier in them, in gen- 
eral. Note, however, that the activated passing over the 
barrier is only one example of kinetically controlled reac- 
tion processes. Another example is the reactions limited 
by jump-like migration with small jump rates at short 
distances (naturally, one should take into account the 
discreteness of space at short distances). In such a case, 
these jumps of small rate can control the rate of reaction. 
The small value of jump rates at short distances r ~ d 
can result from large difference in energy (energy gap) 
between the states of migrating electron and the final 
electron state in the ion pair. Moreover, the rate of jumps 
into this final reacting state (this jumps are probably ir- 
reversible) can be small both for positive and negative 
energy difference (of large absolute value) ^ In the case 
of small values of rates of irreversible jumps into the final 
state the kinetics of the processes is correctly described 
by the Smoluchowski-like equation with partially reflec- 
tive boundary condition, in which the reflective term is 
determined the above-mentioned small (reactive) rate of 
jump in the final state. Analysis of this equation shows 
that the two-state representation of this equation (for the 
Laplace transform R) is also of the form (|2.25p and p.4p 
in the cases of normal diffusion and subdiffusion, respec- 
tively. It is important to emphasize that in the subdif- 
fusion variant of the two-state model (|3.4p the term, de- 
scribing the kinetically controlled reaction in the well, is 
represented in the form e^~"W^", i.e. contains the factor 
e^^" describing the long time memory effects caused by 
anomalous migration. In other words, the kinetic equa- 
tions \2. 25]) and \3.4^ derived in the activated rate model 
appear to he valid for a number of models of kinetically 
controlled reactions in the well, some of which are quite 
applicable to solid state recombination reactions. 

To apply these models one needs to specify the kinetic 
scheme of the process under study which allows one to 
relate the calculated kinetics functions with the observed 
dependences. The scheme was, actually, implied in the 
analysis made in refs. [21] and [26]. We assume that the 
e-h recombination results in the formation of the fluo- 
rescing product (P) whose population will be hereafter 
denoted as np(t). The luminescence intensity I{t) can 
easily be determined from the simple kinetic scheme: 

e + h'-^P ""^^ Po (4.2) 

which takes into account formation of the product P with 
the rate Jr{t) within the well, radiationless deactivation 




FIG. 2: The dependence of SDAR flux Jr{r) (in arbitrary 
units) on dimensionless time t = Wg t for two models of re- 
activity in the well: the activated rate model (full lines) and 
the first order reaction models (squares and triangles). The 
flux is evaluated for: (a) a = 0.3 and (b) a — 0.8, as well as 
for negligibly weak reactivity in the activated rate model and 
for two values of reactivity in the first order reaction model: 
£/ = Wrj,/we^ — 0.001 (trlaugles) and e/ = 0.1 (squares). In 
evaluation the fixed value ja = "ff = 0.2 is used. 

into the ground state Pq with the rate Wd and fluores- 
cence with the rate Wi. In a quite realistic limit of fast 
fluorescence and deactivation when Wd + Wi ^ We, Wr the 
kinetic equation for the population np{t) 

np = ~{wd + 'Wi)np + Jr{t), (4.3) 

corresponding to the scheme (|4.2p . predicts 

I{t) - Winp{t) - [wi/{wd + vu,)]Jr{t). (4.4) 

This formula shows the time dependence of the observed 
luminescence intensity I{t) is proportional to that of the 
recombination (DAR) flux discussed above and therefore 
gives the direct information on the kinetics of geminate 
e-h recombination. 

Noteworthy is that usually in the interpretation of ex- 
perimental results radiationless deactivation is not taken 
into consideration. We mentioned this process only for 
the sake of completeness and because typically it is faster 
than fluorescence transitions. It is seen from the pro- 
posed kinetic scheme that the considered experiments 



10 




10' 10'' 10"' lO"* 10"^ 10'^;(s) 



FIG. 3: Comparison of the time (t) dependence of SDAR flux 
Jr{t) measured experimentally^ (black squares) and calcu- 
lated with eq. (|3.18p in the activated rate (or kinetically con- 
trolled reaction) model of reactivity (full lines) assuming weak 
reactivity, when We^ ~ '"'e„, and for: (1) T — lOS-R', Wg = 
3.3-10*s-\7<, = 0.2; (2) T= 145K, = 1.3 ■ 10^ s'^ ,ya = 
0.3; (3) T = 170A:, w^^ = 3.3 ■ IO^s-^.tI = 0.4. 



can not really give any information about relative effi- 
ciency of these two processes and in what follows we are 
not going to discuss this point. 

Figure 2 demonstrates characteristic features of Jr (t)- 
dependences obtained within the two models of reactiv- 
ity in the well discussed above. It is seen that at in- 
termediate times the first order reaction model predicts 
more sharp changing Jr (t) than the activated rate model, 
though the asymptotic long-time behavior of Jr(t) is the 
same in both models in accordance with derived analyt- 
ical formulas. 

The observed fairly substantial difference between pre- 
dictions of the two models allows for selecting the most 
realistic model by comparison with experimental results. 
Close inspection shows that the very smooth experimen- 
tal dependences Jr {t) are much closer to those predicted 
by the activated rate model rather then by the first order 
reaction one. This means that the activated rate model 
seems to be closer to reality as applied to the process 
under study. 

The fitting of theoretical Jr(t)-dependences to exper- 
imental data (see Fig. 3) enables one to obtain corre- 
sponding parameters of the activated rate model: a, wq^ 
and 7a. Unfortunately, because of the very smooth shape 
of the kinetic dependences these parameters can be de- 
termined fairly approximately. The results can be sum- 
marized as follows: 

(1) The obtained values of a increase approximately 
linearly with the increase of T in agreement with results 



of earlier analysis As for absolute values of the param- 
eter a they are about 10 — 15% lower than those known 
for the same system from time-of-flight experimentsi^ 
This deviation, however, is quite within the accuracy of 
our method (estimated to be ~ 20%). 

(2) The rate wq, (T) increases with temperature T, as 
expected: «7o,(105K) w 3.3 • 10^ s-\ wo^UbK) w 1.3 • 
10^ s-\ and Wo, (170K) « 3.3 • 10^ s-\ 

(3) The parameter 7^ is found to depend on tem- 
perature and ^a{T) increases with T: 7a(105) « 0.2, 
7a(145) « 0.3, and 7a(170) « 0.4 (the estimated ac- 
curacy of these values is about 30%). This behav- 
ior of ^a{T) can be qualitatively interpreted within the 
proposed model. Under quite reasonable assumption 
of small reaction rate in the well, <C We, when 
wg^ « We, , one gets for the ratio 

la(T)ha{To) « {To/T)[weST)/w,STa)r- (4.5) 

In deriving this relation we took into account that 
le{T) /le{To) = Tq/T and assumed that Dg depends on 
temperature weaker than other (rate) parameters. In any 
case the temperature dependence 7a (T) is mainly deter- 
mined by that of We, (T) . 

Estimation [using eq. ((45)) ] yields: 7a(170)/7a(105) « 
1.40 and 7a(145)/7a(105) « 1.31. The corresponding 
values, obtained from fitting kinetic curves, are given 
by: 7a(170)/7a(105) = 1.6 ±0.2 and 7a(145)/7a(105) = 
1.4 ± 0.2. Unfortunately low accuracy of extracted pa- 
rameters of the model does not allow unambiguous test 
of predicted characteristic dependences. 

Concluding this section it is worth discussing the re- 
sults of our analysis and comparing them with those ob- 
tained in the early work [24] (see also ref. [23]), in which 
the kinetics of the same process, geminate e-h recombina- 
tion in amorphous semiconductors a-Si:H, is investigated 
both experimentally and theoretically in a fairly wide re- 
gion of temperatures: 8 < T < 150 K. 

One of principle differences of the interpretation, pre- 
sented this work, from our study lies in the applied model 
of migration of recombining charges in the disordered 
semiconductors. In the work [23] the analysis has been 
made within the Markovian model implying the conven- 
tional Smoluchowski equation of type of eq. (j2.14p for 
the PDF, whereas in our study the migration is assumed 
to be non-Markovian and is described with the use of the 
subdiffusion model. The difference between these two 
descriptions, evidently, shows itself in the long time be- 
havior of the recombination kinetics. Both models pre- 
dict the inverse power type dependence of the reaction 
yield on time Jr{t) ~ 1/^^+"/^, however in the Marko- 
vian case a = 1 while in the non-Markovian (subdiffu- 
sion) one a < 1. Analysis of recent experimental mea- 
surements of the recombination kinetics^^ shows that as 
applied to charge migration in amorphous a-Si:H semi- 
conductors the subdiffusion model is more realistic than 
the Markovian model. This conclusion is strongly sup- 
ported by recent time-of-flight measurements of transient 
currents in these semiconductors/^! Moreover, it is worth 



11 



mentioning that close inspection of experimental results 
presented in the work [24] indicates some deviation of 
the long time part of the recombination kinetics from 
the predicted dependence l/t^^'^. It is also seen that at 
temperatures T > 60 K the long time part is better de- 
scribed by the function l/t^+"/^ with a < 1. 

The intermediate part of the luminescence decay (i.e. 
e — h recombination) kinetics is described in ref. [24] 
by equation of type of eq. ()2.14|) with u{r) = le/r and 
kr ~ e"'"/'"". This complicated equation is solved approx- 
imately by reducing it to the free diffusion one (which in 
turn is solved within the prescribed diffusion approxi- 
mation) and then evaluating the effect of the Coulomb 
potential perturbatively. The weakness of the effect of 
the interaction resulted from the assumed fairly strong 
reactivity at short distances, i.e. fast reactive disappear- 
ance of pairs which is expected to be insignificantly af- 
fected by the potential. This assumption is, however, 
very restrictive because in the process under study the 
Coulomb interaction is quite strong: even at highest 
temperatures studied, T w 150K, for realistic values 
of the contact distance d « lOA and dielectric con- 
stant e « 10 the dimensionless well depth is fairly large, 
u{d) = U{d)l{kBT) « 10. As for lower T « lOK, the 
dimensionless well depth at these temperatures is even 
much (15 times) larger than that at T « 150K. 

The observed reasonable accuracy of the perturba- 
tive solution of the Smoluchowski equation is due as- 
sumed fairly strong and long distant tunneling rate ke ~ 
exp(— r/re), with the tunneling length Ve ~ llA. So 
large tunneling length corresponds to a very small elec- 
tron localization energy Ei « 0.034 eV in amorphous 
a-Si:H, which can easily be estimated by the relation 
re « TB^/EH/Ei^where re « 0.5 A and « 13.5 eV 
are the Bohr radius and the ionization energy of hydrogen 
atom, respectively. The localization energy obtained ap- 
pears to be very close to the thermal energy Eth = 0.026 
eV corresponding to room temperature T = 300 K. Such 
a small value of Ei does not look quite realistic though 
some arguments in favor of this estimate have been pre- 
sented in ref. [24]. 

In our interpretation, instead of assuming very small 
value of El we properly described the effect of the well 
resulting from the attractive Coulomb interaction and 
treated the reactivity suggesting it to be weak enough to 
neglect its manifestation during the time of population 
relaxation within the well. Because of long life time of 
pairs in the well, however, this relatively weak reactiv- 
ity strongly manifests itself in the recombination kinetics 
leading to more smooth recombination kinetics at inter- 
mediate times than that predicted by free diffusion model 
in agreement with experimental results. 

The analysis and comparison with results of experi- 
ments at relatively high temperatures T > 100 K demon- 
strate that the proposed subdiffusion variant of the two- 
state model makes it possible to reproduce the behavior 
of the experimental luminescence decay kinetics (i.e. e- 
h recombination kinetics) both at intermediate and long 



times thus describing fairly accurately the kinetics in the 
whole region of experimentally investigated times. 

Of certain interest is the observed luminescence decay 
kinetics at very low temperature T = 8 K,^^ the slope of 
which in a semilogarithmic scale [i.e. the slope of the de- 
pendence ln/(i)] appeared to change non-monotonically 
with time at intermediate times. In ref. [24] such a be- 
havior is described assuming the recombining pairs to be 
nearly immobile. In this case the kinetics is determined 
by the specific features of kr = kr{r) dependence and 
the initial DPF pi{r) — p{r,t — 0), which in ref. [24] is 
chosen to accurately fit the experimental kinetics. 

It is worth noting that in the proposed two-state model 
similar non-monotonic behavior of the slope of hil(t)- 
dependence is predicted in the first order reaction model 
(see Fig. 2). Such a behavior found without special as- 
sumptions on the initial PDF. However one has to as- 
sume that at low temperatures the first order reaction 
mechanism of reactivity becomes more efficient than the 
kinetically controlled reaction one (which seems to be of 
higher efficiency at high temperatures). In the interpre- 
tation within the first order reaction model one should 
also take into account possible distribution of coordinates 
of the bottom of the well which, for sure, will lead to the 
distribution of average recombination rates. This distri- 
bution can somewhat modify the kinetics as in the case 
of immobile reacting pairs. 



V. CONCLUSIONS 

This work concerns the analysis of the kinetics of gem- 
inate SDARs of interacting particles. The interaction po- 
tential is assumed to be attractive and of the well shape. 
The effect of this interaction is shown to reduce to the 
formation of the quasistatic state within the well (cage) , 
which essentially controls the SDAR kinetics. The reac- 
tion is suggested to occur only in the well. Two models 
of reactivity in the well are discussed: the activated rate 
(or kinetically controlled reaction) and the first order re- 
action models. The results obtained have been used for 
the analysis of the geminate electron-hole recombination 
in amorphous semiconductors a-Si:II. This analysis have 
shown that the first (activated rate) model is able to 
describe the experimental recombination kinetics better 
than the second one and therefore can be considered as 
more appropriate as applied to the processes in amor- 
phous a-Si:H semiconductors at not very low tempera- 
tures. 

Comparison of experimental data with theoretical re- 
sults have enabled us to obtain the values of character- 
istic parameters of the system under study within the 
proposed model. Unfortunately very smooth experimen- 
tal kinetic time dependences does not allow for accurate 
enough determination of the parameters. In addition the 
experimental results presented in refs. [19] and [24] cover 
only relatively narrow temperature region in which pa- 
rameters of the model do not change strongly. 
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For the same reasons it is, strictly speaking, hardly 
possible to make absolutely unambiguous conclusions in 
favor of any of two above-mentioned reactivity models. 
The specific features of the recombination kinetics at 
intermediate times predicted in the first order reaction 
model, which seemingly disagree with those found exper- 
imentally, can, nevertheless, be strongly smoothed in the 
presence of the distribution of kinetic parameters of the 
two-state model (average reaction rate iVr^ and the es- 
caping rate ) . Such a distribution is quite natural for 
studied amorphous semiconductor. 



The results of comparison show that more extensive 
experimental investigations in a wider temperature range 
arc desirable. These studies would provide us with more 
detailed information on the mechanisms of migration and 
interaction of quasiparticles in amorphous semiconduc- 
tors which would allow for making more reliable conclu- 
sions about these mechanisms. 
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